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Abstract. Recent advances in general relativistic magnetohydrodynamic modeling of jets offer 
unprecedented insights into the inner workings of accreting black holes that power the jets in 
active galactic nuclei (AGN) and other accretion systems. I will present the results of recent 
studies that determine spin-dependence of jet power and discuss the implications for the AGN 
radio loud/quiet dichotomy and recent observations of high jet power in a number of AGN. 



1. Introduction 

Relativistic jets are one of the most spectacular manifestations of black hole (BH) accretion. Jet- 
producing accretion systems span 9 decades in central BH mass: from stellar-mass BHs in black 
hole binaries (BHBs) and gamma-ray bursts (GRBs) to supermassive BHs in active galactic 
nuclei (AGN). If a single physical mechanism is responsible for producing jets throughout the 
BH mass spectrum, it must be robust and scale invariant. Magnetic fields are a promising 
agent for jet production because they are abundant in astrophysical plasmas and because the 
properties of magnetically-powered jets scale trivially with BH mass [1-5]. 

How are jets magnetically launched? Figure 1 shows a cartoon depiction of this. Consider a 
vertical magnetic field line attached on one end to a perfectly conducting sphere, which represents 
the central compact object, and on the other end to a perfectly conducting "ceiling" which 
represents the ambient medium (panel a). As the sphere is spinning, after N turns the initially 
vertical field line develops N toroidal field loops (panel b). This magnetic spring pushes against 
the ceiling due to the pressure of the toroidal field. As more toroidal loops form and the toroidal 
field becomes stronger, the spring pushes away the ceiling and accelerates any plasma attached 
to it along the rotation axis, forming a jet (panels (c) and (d) in Figure 1, see the caption for 
details). In the case when the central body is a black hole, which does not have a surface, the 
rotation of space-time causes the rotation of the field lines, and jets form in a similar fashion 
via a process referred to as Blandford-Znajek mechanism (BZ, hereafter) [1]. 




Figure 1. Illustration of jet formation by magnetic fields, (a) Consider a purely poloidal 
(i.e., By = 0) field line attached on one end to a stationary "ceiling" (which represents the 
ambient medium and is shown with hashed horizontal line) and on the other end to a perfectly 
conducting sphere (which represents the central BH or disk and is shown with grey filled circle) 
rotating at an angular frequency Q. (b) After ./V rotations, at time t = t±, the initially purely 
poloidal field line develops iV toroidal loops. This magnetic spring pushes against the "ceiling" 
with an effective pressure p m ~ B^/8tt due to the toroidal field, B v . As time goes on, more 
toroidal loops form, and the toroidal field becomes stronger, (c) At some later time, t = t2, the 
pressure becomes so large that the spring pushes away the "ceiling" and accelerates the plasma 
attached to it along the rotation axis, forming a jet. Asymptotically far from the center, the 
toroidal field is the dominant field component and determines the dynamics of the jet. (d) It is 
convenient to think of the jet as a collection of toroidal field loops that slide down the poloidal 
field lines and accelerate along the jet under the action of their own pressure gradient and hoop 
stress. The rotation of the sphere continuously twists the poloidal field into new toroidal loops 
at a rate that, in steady state, balances the rate at which the loops move downstream. 

2. Radio-Loud/Quiet Dichotomy of Active Galactic Nuclei 

Despite decades of research, we still do not understand what determines jet power in accretion 
systems or how to accurately infer jet powers observationally [6-11]. There are a number of 
unresolved puzzles. For the same nuclear B-band luminosity, radio-loud (RL) AGN have a 
factor of 10 3 higher total (core plus extended) radio luminosity than radio-quiet (RQ) AGN. 
Figure 3 in [12] clearly shows that these two flavors of AGN follow two well-separated tracks 
on the radio-loudness diagram, which shows the radio-loudness parameter, 1Z, the ratio of radio 
to optical luminosities at specified frequencies vs. the Eddington ratio, A. This dichotomy is 
especially pronounced at low A < 0.01 when the accretion systems are thought to be be in a 
radiatively inefficient state of accretion, with a geometrically thick accretion disk. For higher 
A- values, accretion systems can undergo spectral state transitions (as seen in BHBs, [13]), which 
might cause differences in 1Z. Therefore, we limit ourselves to low-A systems. 

If total (core plus extended) radio luminosity is a tracer of jet power, could the dichotomy 
be due to the differences in the spin of central BHs that power the relativistic jets in RL and 
RQ AGN [14]? Characteristic values of BH spin in RQ AGN, which lie predominantly at the 
centers of spiral galaxies, can be quite low, \a\ < 0.3. This can be understood if the orientation 
of angular momentum accreted by the central BHs changes randomly between accretion events 
[15]. On the other hand, the spin of central BHs in RL AGN, which are hosted predominantly 
by elliptical galaxies, can be much higher, a ~ 1 [15]. Can this difference in BH spin lead to a 
factor of 10 3 dichotomy? 

In the presence of large-scale magnetic fields, rotating BHs produce outflows via the BZ 
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mechanism, which extracts BH rotational energy cit ci rsite, 

Pbz = ~, — ^bh o (standard BZ formula, low-spin limit, a 2 C 1), (1) 

where k ~ 0.05 weakly depends on magnetic field geometry, <3?bh is the magnetic flux through 
the BH horizon, and r g = GM/c 2 is BH gravitational radius [1]. This low-spin approximation, 
which we refer to as the standard BZ formula, remains accurate up to a < 0.5 [16, 17], which 
can be seen in Figure 2. Clearly, the oc a 2 scaling does not appear steep enough to explain the 
10 3 dichotomy: BH power varies by a factor of ~ 10 if a varies from 0.3 to 1. Can the power 
dependence becomes steeper as a — > 1? Figure 2 illustrates that an expansion in the powers of 
BH angular frequency, fin = «c/2rn, 

Pbz = — ^^rh^h /(^h) (BZ6 formula, accurate for all values of a), (2) 

remains accurate up to a < 0.95 for / = 1 (referred to as the BZ2 formula), where rn = r g [l + 
(1 - a 2 ) 1 / 2 ] is BH horizon radius. Equation (2) with /(0 H ) « 1 + 1.38{n H r g /c) 2 - 9.2(n H r g /c) A 
remains accurate for all spins [18] , and this gives us our BZ6 formula. Figure 2 shows that while 
at small spin the BZ6 formula (2) and the standard BZ formula (1) agree, as a — > 1 the BZ6 
formula gives about 3 times more power. However, the spin dependence is still not steep enough 
to explain the dichotomy: BH power varies by a factor of ~ 30 if a varies from 0.3 to 1. 

So far we considered razor-thin accretion disks, which is an implicit assumption behind the 
above expressions for power, eqs. (l)-(2). However, disks in low-luminosity AGN are thick, 
with disk angular thickness as large as H/R ~ 1 (including disk body and magnetized corona). 
Such a thick disk mass-loads equatorial field lines on the BH, so that they become part of a 
sub-relativistic disk wind, and only the magnetic flux in the highly magnetized funnel above and 
below the disk contributes to the jet. If the total magnetic flux, $bh, through the BH is held 
constant, this disk occultation effect can cause a much steeper spin-dependence of jet power, 
Pjet oc f2 H , as is seen in Fig. 2. Then, variation of power by a factor of 10 3 is possible, which 
can explain the radio loud/quiet dichotomy of AGN [18]. 

3. What Sets Jet Power in Black Hole Accretion Systems? 

We have shown that BH power is directly proportional to the square of BH magnetic flux, 
$bh, and the square of BH angular frequency, Oh (see eq. 2), with small corrections beyond 
a > 0.95. In Nature, 3>bh is a free parameter, whose value is poorly observationally constrained. 
Clearly, we have ^bh ^ Mc 2 . But what sets the dimensionless ratio, </>g H = ^gjj/Mr^c, which 
characterizes the degree of inner disk magnetization and controls energy extraction from the BH 
[19-24]? Using ^bHj we define BZ efficiency as BZ6 power (eq. 2) normalized by Mc 2 : 

VBZ = ^ x 100% = ^ 2 H ) (^) 2 /(0 H ) x 100%, (3) 

where (...) is a time-average. Previous GRMHD simulations found t?bz ^ 20%, even for nearly 
maximally spinning BHs [25-28]. Are larger values of ??bz possible? This is an especially 
important question since observations suggest that some AGN produce jets and winds at a high 
efficiency, r] = (Pj et + P„i n d) / (Mc 2 ) x 100% > 100% [6-10], where Pj e t and P w i n d are jet and wind 
powers, respectively. Are BHs capable of powering such highly efficient outflows? We tested 
this with global time-dependent GRMHD accretion disk-jet simulations for different values of 
BH spin. As is standard, we initialized the simulations with an equilibrium hydrodynamic torus 
around a spinning BH and inserted a weak, purely poloidal (B^ = 0) magnetic field loop into the 
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Figure 2. Comparison of various approximations for jet power, Pj e t, versus BH angular 
frequency, Oh (lower x-axis), and BH spin, a (upper x-axis). All powers are normalized to the 
maximum achievable power in the standard low-spin BZ approximation, Pj e t oc a 2 , which is 
shown with dot-dashed green line. The standard BZ approximation remains accurate only for 
moderate values of spin, a < 0.5, and for maximally-spinning BHs (a = 1) it under-predicts the 
true jet power by a factor of ~ 3. This can be seen by comparing to the BZ6 approximation, 
Pj e t oc S7 H [l + a(Q}{rg/c) 2 + ^(Oh^/c) 4 ] , which is shown with red solid line and which is 
uniformly accurate for all values of a (the constant factors, a ~ 1.38 and /3 ~ —9.2, are obtained 
in [18]). BZ2 approximation, Pj et oc fi H , shown with blue short dashed line, is accurate up 
to a < 0.95, beyond which it requires a modest correction. If a BH is surrounded by a thick 
accretion disk, the disk intercepts BH power that is emitted into a meridional band occulted by 
the disk. This reduces jet footprint and power relative to the full BH power for razor-thin disks 
given above and shown with red solid line (see also [18]), provided that BH magnetic flux is held 
constant. The resulting BZ6 power for thick disks, with characteristic angular thickness H/R = 1 
expected in low luminosity AGN, is shown with long dashed black curve. Disk occultation effect 
on jet power is most pronounced at low spins: it causes a substantial steepening of jet power 
dependence on spin and can explain a factor of 10 3 radio loud/quiet dichotomy of AGN. 



torus. Clearly, jet efficiency (eq. 3) depends on the magnetic flux, and time-dependent numerical 
simulations show that the larger the large-scale vertical magnetic flux in the initial torus, the 
more efficient the jets [25, 29-31]. To maximize jet efficiency, we populated the torus with a 
much larger magnetic flux than in previous work. In fact, our torus contained more magnetic 
flux than the inner disk can push into the BH. The outcome for BH spin a = 0.99 is shown in 
Figure 3. The magnetorotational instability (MRI, [32]) brings gas and magnetic flux to the 
BH, and an accretion disk of angular density scale height, h/r ~ 0.3, forms. Figure 3(d), (e) 
shows that both (/>bh and rj increase until t ~ 6000r g /c, beyond which they saturate and oscillate 
around the mean, with time-average r\ ~ 140%. At this time, the BH is saturated with magnetic 
flux, and the magnetic field on the hole is so strong that it obstructs the accretion and leads 
to a magnetically-arrested disk, MAD [19, 23, 24, 33-37]. In this state both BH dimensionless 
magnetic flux, 0bh> and outflow efficiency, r], are maximum, so it is not surprising that we find 



4 




* [r„/c] 



Figure 3. A snapshot and time-dependence of a magnetically-arrested accretion disk [19]. 
A movie is available at http://youtu.be/nRGCNaWST5Q . The top panels, (a) and (b), show 
the vertical (x — z) and horizontal (x — y) slices through the accretion flow. Color shows the 
logarithm of density: red shows high and blue shows low values (see color bar). The black circle 
shows a BH (a = 0.99) and black lines show field lines. The bottom 3 panels show, from top 
to bottom, mass accretion rate, M (panel c), dimensionless flux through the BH, (panel 
d), and outflow efficiencies of the whole outflow, r\ (red solid line), the jet, rjj et (green dashed 
line), and the wind, r] w - m d (blue dash-dotted line) (panel e). Since rj ]et and r? w i n d are measured 
at r = 100r 9 , they lag by At ~ W0r g /c relative to rj, which is measured at r = rn- Colored 
symbols show values at the snapshot time, t ~ 12305r 9 /c. At t = 0, the accretion flow contains 
a large amount of large-scale magnetic flux. Accretion brings mass and flux to the hole, and </>bh 
increases until it reaches a maximum value at around t ~ 6000 time units. At this time the BH 
is saturated with magnetic flux and produces as much power as possible. However, the accretion 
flow brings in even more flux, which impedes the accretion and leads to a magnetically arrested 
disk. Some of the flux escapes from the BH via magnetic interchange and flux eruptions, two of 
which are seen in panels (a) and (b), which frees up room for new flux. This process continues 
in a quasi-periodic fashion, and on average the BH produces outflows at rj ~ 140% efficiency, 
i.e., outflows carry more energy than the entire rest-mass energy supplied by accretion. 

much higher values of rj than previously reported [19]. In fact, r/ > 100%, therefore jets and winds 
carry more energy than the entire rest-mass supplied by the accretion, and this unambiguously 
shows that net energy is extracted from the accreting BH. This is the first demonstration of net 
energy extraction from a BH in a realistic astrophysical scenario. Thicker MADs (h/r ~ 0.6) 
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Figure 4. Spin-dependence of various quantities for MADs with h/r ~ 0.3. (a) Dimensionless 
BH spin, 0bh: simulation results (red dots), by-eye fit, 0fi tl comprised of two linear segments 
in a (/>bh _ ^h plane (black line), and 5% uncertainty on the fit (blue band), (b) Energy 
outflow efficiency, rj: simulation results (red dots), BZ6 efficiency (black line, eq. 2) assuming 
</ ) bh(«) = 0fit(a), and 10% uncertainty (blue band), (c) Jet efficiency (green squares) and wind 
efficiency (inverted blue triangles). Dashed line shows 85% of the above BZ6 efficiency (a good 
estimate of jet power for prograde BHs [23]), and blue color shows a 10% uncertainty band, (d) 
BH spin- up parameter, s, for a thin Novikov-Thorne disk [38] (green dash-dotted line) and for 
the simulations (red dots). Whereas for thin disks the equilibrium value of BH spin is a^ T = 1, 
for our simulations it is much lower, is a^ 111 ~ 0.07 (vertical red band). 



produce outflows at an even higher efficiency, r/ ~ 300% [24]. 

Importantly, in the MAD state rj is independent of the initial amount of magnetic flux in 
the accretion flow, i.e., r/ depends only on BH spin, a, and disk density angular thickness, h/r 
[23]. This allows us to reliably study spin-dependence of various quantities, shown in Figure 4. 
Dimensionless BH flux, ^bh ; shows ~ 50% variation with spin (see panel a). Prograde BHs 
(a > 0) are a few times more efficient than retrograde BHs (a < 0) for the same value of |a| 
(see panel b). Quite the opposite — 10 times larger efficiency of retrograde BHs than of prograde 
BHs — is predicted by the "gap" model [39], which assumes that zero magnetic flux populates 
the "gap" between the BH horizon and the innermost stable circular orbit (ISCO). We find that 
this interpretation of the ISCO — as the surface inside of which all magnetic flux nearly freely 
falls into the BH — does not hold in MADs: interchange instability [40] and flux eruptions (see 
Fig. 3) populate the ISCO interior with flux. In fact, the thinner the disk, the larger the fraction 
of BH magnetic flux that resides in the region between the BH horizon and the ISCO [23, 24]. 

Panel (c) shows the division of total outflow efficiency into highly magnetized jet and weakly 
magnetized wind components, with efficiencies r/j et and r/ w i n d, respectively. Since jets are BH 
spin-powered (eq. 2), for a = jet efficiency vanishes, but winds still derive their power from 
an accretion disk via a Blandford-Payne-type mechanism [41]. The larger the spin, the more 
efficient jets and winds. For rapidly spinning BHs most of the energy — about 85% for prograde 
BHs — is carried by relativistic jets. Importantly, even rather slowly spinning BHs, with a < 0.5, 
produce prominent BH spin-powered jets. This is in agreement with the recent evidence that 



6 



jets in BHBs are powered by BH rotation over a wide range of BH spin, 0.1 < a < 0.9 [42]. 

Do our highly efficient jets affect the spin of central BHs? Figure 4(d) shows spin-dependence 
of BH spin-up parameter, s = M/M x da/dt [43]. For standard geometrically thin accretion 
disks, s > (see Figure 4(d) and [38]), and the equilibrium spin is a[I T = 1 [44] (we neglect 
photon capture by the BH, which would limit the spin to a ~ 0.998 [45]). Thick accretion flows 
in time-dependent numerical simulations [29, 43, 46] and semi- analytic studies [47] typically 
have o eq ~ 0.9 (however, see [48]). Figure 4(d) shows that the equilibrium value of spin for our 
MADs (with h/r ~ 0.3) is much smaller, ~ 0.07, due to large BH spin-down torques by 
our powerful jets. Since jet efficiency is ~ 10 3 times lower at this value of spin than for rapidly 
spinning BHs, such equilibrium spin systems are possible candidates for RQ AGN. 

4. Conclusions 

We confirm that the standard BZ power formula, Pj e t oc a 2 (eq. 1), remains accurate only for 
a < 0.5 and present a higher order BZ6 approximation (eq. 2) that is accurate for all values of 
spin (see Figure 2). We find that if the BH magnetic flux is held constant, the presence of a 
thick accretion disk in low-luminosity AGN leads to a steep dependence of BH jet power on spin, 
.Fjet oc fi^j. This steep dependence allows to explain a factor of 10 3 radio loud/quiet dichotomy 
of AGN via having two galaxy populations different by the spin of central BHs. 

We carried out a series of time-dependent global GRMHD simulations of BH accretion that 
contain a large amount of magnetic flux. We show that the accumulation of magnetic flux 
around the center can lead to the formation of a strong centrally-concentrated magnetic field 
that saturates the BH, obstructs the accretion, and leads to a magnetically- arrested disk (MAD, 
[19, 23, 24, 33-37]). We show that in this state the outflow efficiency, rj, depends only on the BH 
spin, a, and the angular density thickness of the accretion flow, h/r, and is independent of the 
initial amount of magnetic flux in the disk. Since BH magnetic flux is as large as possible, we 
expect MAD systems to achieve the maximum r/ for a given BH spin and disk thickness. Indeed, 
we find highly efficient outflows, with rj > 100%, which suggests that MADs might explain 
observations of AGN with apparent r\ ~ few x 100% [6-11]. We determine the spin-dependence 
of jet, wind, and total efficiencies in our simulations (Figure 4) and expect that this information 
can be useful for calibration of semi-analytic jet power models, for studies seeking to infer BH 
spin from the observed jet efficiency, and for tests of general relativity [11, 49-52]. 
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